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ABSTRACT 

Change of sign of the LNRF-velocity gradient has been found for accretion 
discs orbiting rapidly rotating Kerr black holes with spin a > 0.9953 for 
Keplerian discs and a > 0.99979 for marginally stable thick discs. Such 
a "humpy" LNRF-velocity profiles occur just above the marginally stable 
circular geodesic and could be related to oscillations of accretion discs. The 
frequency of such "hump" -induced oscillations can be identified with the 
maximal rate of change of the orbital velocity within the "humpy" profile. 
Therefore, we introduce an extended orbital resonance model (EXORM) of 
quasiperiodic oscillations (QPOs) assuming non-linear resonant phenomena 
between oscillations with the orbital epicyclic frequencies and the humpy 
frequency defined in a fully general relativistic way. The EXORM is de- 
veloped for both Keplerian discs and perfect-fluid tori where the approxi- 
mation of oscillations with epicyclic frequencies is acceptable. Clearly, the 
EXORM could be applied to the near-extreme Kerr black hole systems ex- 
hibiting relatively complex QPO frequency patterns. Assuming a Keplerian 
disc, it can be shown that in the framework of the EXORM, all the QPOs 
observed in the microquasar GRS 1915+105 could be explained, while it 
is not possible in the case of QPOs observed in the Galactic Centre source 
SgrA*. 
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1 INTRODUCTION 

High frequency (kHz) twin peak quasi-periodic oscillations (QPOs) with frequency 
ratios 3 : 2 (and sometimes 3:1) are observed in microquasars (see, e.g., van der 
Klis, 2000; McClintock and Remillard, 2004; Remillard, 2005). In the Galactic 
Centre black hole SgrA*, Genzel et al. (2003) measured a clear periodicity of 
1020 sec in variability during a flaring event. This period is in the range of Ke- 
plerian orbital periods at a few gravitational radii from a black hole with mass 
M ~ 3.6 x 10 6 M Q estimated for SgrA* (Ghez et al, 2005; Weinberg et al., 2005). 
More recently, Aschenbach et al. (2004); Aschenbach (2004, 2006) reported three 
QPO periodicities at 692 sec, 1130 sec and 2178 sec that correspond to frequency 
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ratios (1/692) : (1/1130) : (1/2178) - 3:2:1. These observational data are not quite 
convincing (see, e.g., Abramowicz ct al., 2004), but surely deserve attention (As- 
chenbach, 2007). 

Detailed analysis of the variable X-ray black-hole binary system (microquasar) 
GRS 1915+105 reveals high-frequency QPOs appearing at five frequencies, namely 
vi = (27 ± 1) Hz (Belloni et al., 2001), v 2 = (41 ± 1) Hz, v 3 = (67 ± 1) Hz (Morgan 
et al., 1997; Strohmayer, 2001), and v± = (113 ±5) Hz, v 5 = (167 ±5) Hz (Remil- 
lard and McClintock, 2006). In this range of their errors, both upper pairs 
are close to the frequency ratio 3 : 2 suggesting the possible existence of reso- 
nant phenomena in the system. Observations of oscillations with these frequen- 
cies have different qualities, but in all five cases the data are quite convinc- 
ing (see McClintock and Rcmillard, 2004; Rcmillard and McClintock, 2006). 

Several models have been developed to explain the kHz QPO frequencies, and it 
is usually preferred that these oscillations are related to the orbital motion near the 
inner edge of an accretion disc. In particular, two ideas based on the strong-gravity 
properties have been proposed. While Stella and Vietri (1998, 1999) introduced the 
"Relativistic Precession Model" considering that the kHz QPOs directly manifest 
the modes of a slightly perturbed (and therefore epicyclic) relativistic motion of 
blobs in the inner parts of the accretion disc, Kluzniak and Abramowicz (2001) 
propose models based on non-linear oscillations of an accretion disc that assume 
resonant interaction between orbital and/or epicyclic modes. In a different context, 
the possibility of resonant coupling between the epicyclic modes of motion in the 
Kerr spacetime was also mentioned in the early work of Aliev and Galtsov (1981). 
The radial and vertical epicyclic oscillations could be related to both the thin Ke- 
plerian discs (Abramowicz et al., 2003; Kato, 2001) and the thick, toroidal accre- 
tion discs (Rezzolla et al., 2003). In particular, the observations of high frequency 
twin peak QPOs with the 3 : 2 frequency ratio in microquasars can be explained 
by the parametric resonance between the radial and vertical epicyclic oscillations, 
v y : v T — 3:2. This hypothesis, under the assumption of geodesic oscillations (i.e., 
for thin discs), puts strong limit on the mass-spin relation for the central black hole 
in microquasars (Torok et al., 2005; Torok, 2005; Torok et al., 2006). 

Aschenbach (2004, 2006, 2007) discovered that two changes of sign of the radial 
gradient of the Keplerian orbital velocity as measured in the locally non-rotating 
frame (LNRF, Bardeen et al., 1972) occur in the equatorial plane of Kerr black holes 
with a > 0.9953. Stuchlik et al. (2005) have found that the gradient sign change 
in the LNRF-velocity profiles occurs also for non-geodesic motion with uniform 
distribution of the specific angular momentum £(r,6) = const (i.e., in marginally 
stable thick discs) around extremely rapid Kerr black holes with a > 0.99979. 1 The 
global character of the phenomenon is given in terms of topology changes of the von 
Zeipel surfaces (equivalent to equivelocity surfaces in the tori with £(r, 9) = const). 



1 Note that the assumption of uniform distribution of the specific angular momentum can be 
relevant at least at the inner parts of the thick disc and that matter in the disc follows nearly 
geodesic circular orbits nearby the center of the disc and in the vicinity of its inner edge determined 
by the cusp of its critical equipotcntial surface (see Abramowicz et al., 1978). 
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Toroidal von Zeipel surfaces exist around the circle corresponding to the minimum 
of the equatorial LNRF- velocity profile, indicating possibility of development of 
some instabilities in that part of the marginally stable disc with positive gradient 
of the orbital velocity in LNRF (Stuchlfk ct al., 2004, 2005, 2006a,b, 2007a,b,c,d). 

The positive radial gradient of orbital LNRF- velocity around black holes with a > 
0.9953 seems to be a physically interesting phenomenon, even if a direct mechanism 
relating this phenomenon to triggering the oscillations, and subsequent linking of 
the oscillations to the excitation of radial (and vertical) epicyclic oscillations, is 
unknown. Therefore, an extended orbital resonance model (EXORM) has been 
developed, with hypothetical hump-induced oscillations assumed to enter a non- 
linear resonance with the radial or vertical epicyclic oscillations (Stuchlfk et al., 
2007b). It should be stressed that due to the non- linear resonance, combinational 
frequencies are allowed to be observable. 

In the EXORM, the frequency of the hump-induced oscillations is related to the 
maximal positive radial gradient of the LNRF-velocity in the "humpy" velocity 
profile in the general relativistic, coordinate-independent form. Further, since the 
gradient is defined locally, being connected to the LNRF, it has to be transformed 
into the form related to distant stationary observers, giving observationally rele- 
vant "humpy" frequency Uh. Then the "humpy" and epicyclic frequencies could be 
estimated at the radius of definition of the "humpy" frequency. 

In the case of Keplerian discs, the epicyclic resonance radii r^-.i and r^-.x (with 
v v : v x = 3:1, 4:1) are located in vicinity of the "humpy" radius where efficient 
triggering of oscillations with frequencies ~ could be expected. Asymptotically 
(for 1 — a < 10~ 4 ) the ratio of the epicyclic and Keplerian frequencies and the 
humpy frequency is nearly constant, i.e., almost independent of a, being for the 
radial epicyclic frequency v T : <~ 3 : 2. In the case of thick discs, the situation is 
more complex due to dependence on distribution of the specific angular momentum I 
determining the disc properties. For 1 — a < 10~ 6 , the frequency ratios of the humpy 
frequency and the orbital and epicyclic frequencies arc again nearly constant and 
independent of both a and I being for the radial epicyclic frequency v T : z/h ~ 4:1. 
In the limiting case of very slender tori [l <~ £ ms ) the epicyclic resonance radius 
i"4:i ~ rh for all the relevant interval of 1 — a < 2 x 10~ 4 . 

In Section 2, we briefly summarize properties of the Aschenbach effect for Kep- 
lerian thin discs, and I = const thick discs. In Section 3, the extended resonance 
model is introduced, i.e., the critical "humpy" frequency, connected to the LNRF- 
velocity positive gradient in the humpy profiles, is given in the physically relevant, 
coordinate independent form for the both Keplerian and t = const discs. At the 
radius of its definition, the humpy frequency is compared to the radial and vertical 
epicyclic frequency and the orbital frequency. In Section 4, fitting of the observed 
frequencies in the GRS 1915+105 microquasar in the framework of the EXORM is 
summarized, while it is demonstrated that the data reported for Sgr A* could not 
be probably fitted by EXORM. In Section 5, concluding remarks are presented. 
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2 LNRF- VELOCITY PROFILES OF ACCRETION DISCS 

The locally non-rotating frames (LNRF) are given by the tetrad of 1-forms (Bardeen 
et al., 1972) 



: (*) = V2 dt , e (<p) = (J^j V2 sin 9 {Ap -udt). 



(1) 



where 

» = = (2) 

9<p<p A 

is the angular velocity of the LNRF relative to distant observers. 

In the Kerr spacetimes with the rotational parameter assumed to be a > 0, the 
relevant metric coefficients in the standard Boyer-Lindquist coordinates read: 

A-a 2 sin 2 6» 2arsin 2 _ A sin 2 6 _£ 

9tt — ^ , gtip ~ ^ , 9iptp — — ^ — , 9rr — , 9ee — ^ j (3) 

where 

A = r 2 - 2r + a 2 , S = r 2 + a 2 cos 2 9, A = (r 2 + a 2 ) 2 - Aa 2 sin 2 6 . (4) 

The geometrical units, c = G = 1, together with putting the mass of the black hole 
M = 1, are used in order to obtain completely dimensionless formulae hereafter. 

For matter orbiting a Kerr black hole with a 4-velocity and angular velocity 
profile f2(r, 0), the azimuthal component of its 3- velocity in the LNRF reads 

v<*> = ^ = ^2(n- w) . (5) 



2.1 Keplerian thin discs 

In thin discs matter follows nearly circular geodetical orbits characterized by the 
Keplerian distributions of the angular velocity and the specific angular momentum 
(in the equatorial plane, 9 = n/2) 

o o i \ 1 o o i \ r2 - 2arl/2 + fl2 fR \ 

Q = K (r;a)= (r 3 /2+a) , * = l K ( r;o) = . (6) 

The azimuthal component of the Keplerian 3-velocity in the LNRF reads 
(v) _ (r 2 +a 2 ) 2 -a 2 A-2ar(r 3 / 2 + a) 
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and formally diverges for r — ► r + = 1 + \/I — a 2 , where the black-hole event horizon 
is located. Its radial gradient is given by 

dvjf } _ r 5 + a 4 (3r + 2) - 2aV/ 2 (3r + 1) - 2a 2 r 2 (2r - 5) + 2ar 5 / 2 (5r - 9) 
dr ~ 2A 3 / 2 v^(r 3 / 2 + a) 2 

(8) 

As shown by Aschenbach (2004, 2006), the velocity profile has two changes of the 
gradient sign (where dV^ /dr = 0) in the field of rapidly rotating Kerr black holes 
with a > a c (K) = 0.9953. 



2.2 Marginally stable tori 

Perfect-fluid stationary and axisymmetric toroidal discs are characterized by the 
4-velocity field = (U\ 0, 0, Uf) with £/* = XP{r,ff), U v = U v {r,9), and by 
distribution of the specific angular momentum £ = —U v /Ut- The angular velocity 
of orbiting matter, Q = U v /£/*, is then related to £ by the formula 

n = _tgtt + 9t v (9) 
Zgtf + g vv 

The marginally stable tori are characterized by uniform distribution of the specific 
angular momentum 

I = e(r, 0) = const , (10) 

and are fully determined by the spacetime structure through equipotcntial surfaces 
of the potential W — W(r, 9) defined by the relations (Abramowicz et al., 1978) 

(U t )in ' " Qui 2 + 2g tv £ + g vv ' 

the subscript "in" refers to the inner edge of the disc. The LNRF orbital velocity 
of the torus is given by 

V ( V ) = A(A-a 2 sm 2 9) + 4a 2 r 2 sin 2 9 £ 
T SVA(A- 2afr)sin6» 

For marginally stable tori it is enough to consider the motion in the equatorial 
plane, 9 — ir/2. Formally, this velocity vanishes for r — > oo and r — > r+, i.e., there 
must be a change of its radial gradient for any values of the parameters a and £, 
contrary to the case of Keplerian discs. The radial gradient of the equatorial LNRF 
velocity of I = const tori reads 

dV^ ] _ f [A + (r - l)r][r(r 2 + a 2 ) - 2a(£ - a)} 



W ^ W . n = l n J^ 7 {Ut f = y gttg ^ . (n) 



dr \ [r(r 2 + a 2 )-2a(£-a)] 2 VA 



r (3r 2 + a 2 )A \ 



[r(r 2 + a 2 ) - 2a(£ - a)] 2 VA J 



£, (13) 
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so it changes its orientation at radii determined for a given £ by the condition 
1 = 1 (r-a)-al ^ +a a )(r - 1) - 2rA] 

* *«*ir> a )- a + 2a[A + r(r-l)] ' 1 ' 

For both thick tori and Keplerian discs we have to consider the limit on the 
disc extension given by the innermost stable orbit. For Keplerian discs this is 
the marginally stable geodetical orbit, r in r ms , while for thick tori this is an 
unstable circular geodesic kept stable by pressure gradients and located between 
the marginally bound and the marginally stable geodetical orbits, r m b < r- ln < r ms , 
with the radius being determined by the specific angular momentum £ = const £ 
(^ms)^mb) through the equation £ = £x(r;a); C (4nb) denotes specific angular 
momentum of the circular marginally stable (marginally bound) geodesic. 

Detailed discussion of Stuchlik et al. (2005) shows that two physically relevant 
changes of sign of dV^ /dr in the tori occur for Kerr black holes with the ro- 
tational parameter a > a c (T) = 0.99979. The interval of relevant values of the 
specific angular momentum I G (f ms (a),4 x ( m ax)( fl ))i where ^ C x(max)( a ) corresponds 
to the local maximum of the function (14), grows with a growing up to the crit- 
ical value of a c ( mb ) = 0.99998. For a > a c ( m b), the interval of relevant values 
of £ e (4is(a),Cb(a)) is narrowing with the rotational parameter growing up to 
a = 1, which corresponds to a singular case where £ ms (a = 1) = £ m b(a = 1) = 2. 
Notice that the situation becomes to be singular only in terms of the specific angu- 
lar momentum; it is shown (see Bardeen et al., 1972) that for a — 1 both the total 
energy E and the axial angular momentum L differ at r ms and r m b, respectively, 
but their combination, £ = L/E, giving the specific angular momentum, coincides 
at these radii. 

A physically reasonable global quantity characterizing rotating fluid configura- 
tions in terms of the LNRF orbital velocity is so-called von Zcipel radius defined by 
the relation 

ft=-*= (1-^)0, (15) 
''lnr.f 

which generalizes in another way as compared with (Abramowicz et al., 1995) the 
Schwarzschildian definition of the gyration radius g (Abramowicz et al., 1993). Note 
that, except for the Schwarzschild case a = 0, the von Zcipel surfaces, defined as the 
surfaces of lZ(r, 8; a, £) = const, do not coincide with those introduced by Kozlowski 
et al. (1978) as the surfaces of constant £/VL (see Stuchlik ct al., 2005 for more 
details). 

In the case of marginally stable tori the von Zeipel surfaces 1Z = const coincide 
with the equivelocity surfaces V^(r, 0; a, £) = v!f> = const. Topology of the von 
Zeipel surfaces can be directly determined by the behaviour of the von Zeipel radius 
in the equatorial plane 

^, „ /„ „x r(r 2 + a 2 ) - 2a(£- a) ,,„. 

Tl(r,0 = n/2;a,£) = ^ i ~ 7 =-^ '- . (16) 

rv A 



Extended resonance model 455 




0.5 1 1.5 2 2.5 3 

r sin 




r sin 9 



Figure 1. Von Zeipel surfaces (meridional sections). For a > a c ( T ) and £ appropriately 
chosen, two surfaces with a cusp, or one surface with both the cusps, together with closed 
(toroidal) surfaces, exist, being located always inside the ergosphere (dashed surface) of a 
given spacetime. Both the outer cusp and the central ring of closed surfaces are located 
inside the toroidal equilibrium configurations corresponding to marginally stable thick 
discs (light-gray region; its shape is determined by the critical self-crossing equipotential 
surface of the potential W(r, 0). The cross (+) denotes the centre of the torus. Dark 
region corresponds to the black hole. Figures illustrating all possible configurations of the 
von Zeipel surfaces are presented in Stuchlfk et al. (2005). Here we present the figure 
plotted for the parameters a — 0.99998, £ = 2.0065. Critical value of the von Zeipel 
radius corresponding to the inner and the outer self-crossing surface is 72. c ( in ) = 3.429 and 
7?-c(out) = 3.804, respectively, the central ring of toroidal surfaces corresponds to the value 
7^-contcr = 3.817. Interesting region containing both the cusps and the toroidal surfaces is 
plotted in detail at the left lower figure. Right lower figure shows the behaviour of the von 
Zeipel radius in the equatorial plane. (Taken from Stuchlfk et al., 2007b) 
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The local minima of the function (16) determine loci of the cusps of the von Zcipel 
surfaces, while its local maximum (if it exists) determines a circle around which 
closed toroidally shaped von Zeipel surfaces are concentrated (see Fig. 1). Notice 
that the inner cusp is always physically irrelevant being located outside of the 
toroidal configuration of perfect fluid. Behaviour of the von Zeipel surfaces nearby 
the centre and the inner edge of the thick discs orbiting Kerr black holes with 
a > a c(T) = 0.99979, i.e., the existence of the von Zcipel surface with toroidal 
topology, suggests possible generation of instabilities in both the vertical and radial 
direction. 

In terms of the redefined rotational parameter (1 — a), the "humpy" profile of 
the LNRF orbital velocity of marginally stable thick discs occurs for discs orbiting 
Kerr black holes with 1 — a < 1 — a c rr) = 2.1 x 1CP 4 , which is more than one order 
lower than the value 1 — a c (K) = 4.7 x 1CP 3 found by Aschenbach (2004) for the 
Keplerian thin discs. Moreover, in the thick discs, the velocity difference AV^ is 
smaller but comparable with those in the thin discs. In fact, for a — > 1, the velocity 
difference in the thick discs AV^ w 0.02, while for the Keplerian discs it goes even 
up to AV^ } w 0.07 (Stuchlik ct al., 2007d). 



3 EXTENDED ORBITAL RESONANCE MODEL 

The orbital resonance model assumes non-linear parametric and forced resonances 
of oscillations with the orbital (Keplerian) and radial or vertical epicyclic frequen- 
cies, or their combinations. Here, we extend this model by introducing hypothetical 
additional oscillations, induced by the hump in the LNRF-velocity profile of accre- 
tion discs of both Keplerian and toroidal character, that are supposed to be in a 
non- linear resonance with orbital or epicyclic oscillations. 

In Kerr spacetimes, the frequencies of the radial and latitudinal (vertical) epicyclic 
oscillations related to an equatorial Keplerian circular orbit at a given r are deter- 
mined by the formulae (see, e.g., Aliev and Galtsov, 1981) 



where the Keplerian frequency z/k = Ok/27t. A detailed analysis of properties 
of the epicyclic frequencies can be found in Torok and Stuchlik (2005a, b). The 
epicyclic oscillations with the frequencies v T , u v can be related to both the thin 
Keplerian discs (Abramowicz and Kluzniak, 2001; Kato, 2004) and thick, toroidal 
discs (Rezzolla et al., 2003). 

According to Aschenbach (2004, 2006), the non-monotonicity of the LNRF-velocity 
profile of accretion discs could excite oscillations with characteristic frequency that 
has to be related to the maximum gradient in the "humpy" part of the accretion 
discs velocity profile. 

Although there is no detailed idea on the mechanism generating the hump-induced 
oscillations, it is clear that the Aschenbach proposal of defining the characteristic 



v 2 = i4(l - Gr- 1 + 8ar- 3/2 - Za 2 r- 2 ) 
is 2 = is 2 = 4(1 - Aar- 3/2 + 3a 2 r~ 2 ) , 



(17) 
(18) 
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frequency deserves attention. It should be stressed, however, that a detailed analysis 
of the instability could reveal a difference between the characteristic frequency and 
the actual observable one, as the latter should be associated with the fastest growing 
unstable mode. In any case, the humpy frequency represents an upper limit on the 
frequencies of the hump-induced oscillations, as it is given by maximum of the 
LNRF-velocity gradient in the humpy part of the velocity profile. 

At the present state of the EXORM, we assume that the characteristic humpy fre- 
quency is a typical frequency of oscillations induced by the conjectured "humpy in- 
stability," and that the humpy oscillations could excite oscillations with the epicyclic 
frequencies or some combinational frequencies, if appropriate conditions for a forced 
resonance are satisfied in vicinity of the radius where the humpy oscillations oc- 
cur (Stuchhk ct al., 2007d). 

The fully general relativistic definition of the critical frequency for a possible 
excitation of oscillations in the disc is given by the relations 



R _ dVM 

IS--;*. 



OR 




(19) 



where = (r; a) in thin Keplerian discs, and = (r; I, a) in marginally 
stable thick discs and R is the physically relevant, coordinate independent proper 
radial distance. Such a locally defined frequency, confined naturally to the ob- 
servers orbiting the black hole with the LNRF, should be further related to distant 
stationary observers by the formula (taken at the BL coordinate r corresponding to 
(dV^/dR) max ) 



Uh = = yj-(g tt + 2ugt<p + u 2 g vv ) i/£ it . (20) 

We call such a coordinate-independent and, in principle, observable frequency the 
"humpy frequency," as it is related to the humpy profile of V^', and denote it 
v^. It should be stressed that the physically relevant humpy frequency = z/^, 
connected to observations by distant observers and exactly defined by Eqs (19) and 
(20), represents an upper limit on characteristic frequencies of oscillations induced 
by the hump of the LNRF-velocity profile, and the realistic humpy frequencies, 
as observed by distant observers, can be expected close to but smaller than v^. 
Further, we denote the BL radius of definition of the humpy oscillations frequency, 
where dV^/dR = {dV^ /dR) max . 

In the case of the Keplerian discs we obtain the "humpy frequency" to be given 
by the relation 

_ -rj - a 4 (3r h + 2) + 2a 3 r^ /2 (3r h + 1) - 2a 2 r 2 (2r h - 5) + 2ar 5 h /2 (5r h - 9) 

2A h r2 (r 3 h /2 + a) 2 



/ 4a 2 
h_2 ~r h (r 2 +a 2 ) + 2a 2 ' (21) 
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where A h = r\ — 2r h + a 2 . The BL radius r h where the positive gradient of the 
velocity profile in terms of the proper radial distance reaches its maximum, so-called 
"humpy radius," is given by the condition 

leading to the equation 

3a 7 (r + 2) + a 6 y/r{21r 2 + 18r - 4) - a 5 r(33r 2 + lOr + 20) 

+ a 4 ry/r{45r 3 - 62r 2 - 68r + 16) - a 3 r 3 (83r 2 - 122r - 60) 

+ a 2 r 4 V^(27r 2 - 130r + 136) - 9ar 5 (7r 2 - 26r + 24) 

+ r 7 V?(3r - 2) = , (23) 

which must be solved numerically. The spin dependence of the humpy radius and 
the related humpy frequency is illustrated in Fig. 2. The humpy radius rh falls 
monotonically with increasing spin a, while the humpy frequency has a maxi- 
mum for a = 0.9998, where ^ n ( max ) = 607 (M Q /M) Hz, and it tends to ^h(a^i) = 
588 (M©/M) Hz. 

The ratios of the humpy frequency and the orbital and epicyclic frequencies at 
the humpy radius were determined in Stuchlik et al. (2007b) revealing almost spin- 
independent asymptotic behaviour for a — > 1 represented closely by the ratios of 
integer numbers, vk : v w : v Y : ~ 46 : 11 : 3 : 2, which imply a possibility of resonant 
phenomena between the hump-induced and orbital or epicyclic oscillations. For 
Kcplerian discs, the ratios of the epicyclic frequencies and the humpy frequency are 
given in the dependence on the black- hole spin in Fig. 3. 

The marginally stable tori have a structure that depends on the value of the 
specific angular momentum £ G (^ m s7^mb)- The oscillations of slender tori [l w 
4n S ) have frequencies equal to the epicyclic frequencies relevant for test particle 
motion, but the frequencies of non-slender tori are different, as shown for pseudo- 
Newtonian tori (Sramkova, 2005; Blaes et al., 2007) and expected for tori in the 
strong gravitational field of Kerr black holes. Therefore, comparison of the humpy 
frequencies and the epicyclic frequencies is relevant for the slender tori only. 

The humpy frequency is defined for all a > 0.99979 and all t € (^ms,^mb), see 
Fig. 4. It is important that in the field of Kerr black holes with 1 — a < 10~ 8 , there is 
z/h(a, t) ~ 150 Hz (M/Mq)' 1 independently of a and £ (Stuchlik et al., 2007b). Fur- 
ther, the physically important case of tori admitting evolution of toroidal von Zcipel 
surfaces with the critical surface self-crossing in both the inner and the outer cusps 
is allowed at i = € cr i t , where £ cr - lt > £ ms only slightly differs from £ ms , i.e., such tori 
can be slender, see Fig. 4. The ratios of v^jv^, v^jv^ and ^ /^h ar e given for the 
tori with I w £ ms in Fig. 5. Their asymptotical values, valid for 1 — a < 10~ 6 , are 
independent of both a and £. 

Of course, in realistic situations the hump-induced oscillation mechanism could 
work at the vicinity of rh, with slightly different frequencies; we should take into 
account that the shift of the radius, where the mechanism works, shifts both the 
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Figure 2. Spin-dependence of the humpy frequency and the humpy radius rh that is 
compared with the Boyer-Lindquist radius of the innermost stable circular orbit. (Taken 
from Stuchlfk et al., 2007d) 



locally measured (LNRF) frequency (Eq. (19)) and the frequency related to distant 
observers (Eq. (20)). The zones of radii, where the critical frequency v^ it differs 
up to 1 %, 10% and 20% of its maximal value (given by (dV^ / dR) max ) for thin 
(Keplerian) discs or 1 %, 5% and 10% of its maximum for marginally stable discs 
with t = £ ms , are given in Fig. 6. We can see (Fig. 6) that the resonant epicyclic 
frequencies radii r&i and r±.\ are located within the zone of the hump-induced 
oscillation mechanism in both thin discs and marginally stable tori. 

In Keplerian discs the sign changes of the radial gradient of the orbital velocity 
in LNRF occur nearby the r = r 3: i orbit (with v v : v T = 3 : 1), while in the 
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Figure 3. Spin dependence of the ratios of the radial (i/ r ) and vertical (i/ v ) epicyclic 
frequencies, and the Keplerian frequency (uk) to the thin-disc humpy frequency related 
to distant observers {vh). All the frequency ratios are asymptotically (for 1 — a < 1CP 4 ) 
constant. There si uk '■ '■ Vr ■ Vh ~ 46 : 11 : 3 : 2. Therefore, we can expect some resonant 
phenomena on the ratio of u r : Vh ~ 3 : 2, and i^k : fv ~ 4 that could be both correlated. 
(Taken from Stuchh'k et al., 2007b) 
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Figure 4. Interval of humpy frequencies for the marginally stable thick discs with £ £ 
(^ms,^mt>) as a function of the black-hole spin a. For a-tl, the interval is narrowing and 
asymptotically reaching the value of 150 Hz (M/Mq) -1 . Dotted curve corresponds to the 
humpy frequencies of marginally stable slender tori with £ = £ cr it, for which the critical 
von Zeipel surface contains two cusps (as it is demonstrated for one special case in the left 
panel of the figure; the torus is given by the light-gray region). 
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Figure 5. Spin dependence of the ratios of the radial (u T ) and vertical (v v ) epicyclic 
frequencies, and the orbital frequency (v ) of the marginally stable I — £ m s disc to the 
thick-disc humpy frequency related to distant observers (I'h)- All the frequency ratios are 
asymptotically (for 1 — a < 1CP 6 ) almost constant. (Taken from Stuchh'k et al., 2007b) 
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Figure 6. Positions of local extrema of (in BL coordinates) for Keplerian discs (a) 
and marginally stable discs with £ = £ ms (b) together with the locations of resonant orbits 
rs-.i and taa (where the resonance between the vertical and radial epicyclic oscillations 
takes place) in dependence on the rotational parameter a of the black hole. Dashed curve 
corresponds to the maximum positive values of the LNRF orbital velocity gradient in 
terms of the proper radial distance where the critical frequency u^ it is defined, boundaries 
of shaded regions correspond to orbits where the velocity gradient giving the character- 
istic frequency, dV M /dR, reaches (a) 99%, 90%, 80% and (b) 99%, 95%, 90% of its 
maximum. (Taken from Stuchh'k et al., 2007b) 



vicinity of the r = r 3:2 orbit (with v v : v r = 3 : 2), dV^/dr < for all values 
of a for both Keplerian discs and marginally stable tori with all allowed values 
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of I. The parametric resonance, which is the strongest one for the ratio of the 
epicyclic frequencies v w : v T = 3 : 2, can occur at the r = r 3:2 orbit, while its effect 
is much smaller at the radius r = r^i, as noticed by Abramowicz ct al. (2003). 
Nevertheless, the forced resonance may take place at the r^-.i orbit. Notice that the 
forced resonance at r = r 3: i can generally result in observed QPOs frequencies with 
3 : 2 ratio due to the beat frequencies allowed for the forced resonance as shown 
in Abramowicz et al. (2004). But the forced resonance at r^.i between the epicyclic 
frequencies, induced by the humpy profile of V^, seems to be irrelevant in the case 
of microquasars, since all observed frequencies lead to the values of the rotational 
parameter a < a c ( K ), as shown by Torok et al. (2005). 



4 APPLICATION OF THE EXTENDED RESONANCE MODEL 

The extended orbital resonance model with hump-induced oscillations can be ap- 
plied only to black hole systems containing a near-extreme Kerr black hole candi- 
dates. One of the most promising of such systems seems to be the microquasar 
GRS 1915+105, where the extremely high spin a ~ 1 was predicted by continuous 
spectra fitting method (McClintock et al., 2006). Another promising candidate for 
the near-extreme Kerr black hole could be considered the Galaxy Centre SgrA*. 
In fact, all the QPO frequencies observed in GRS 1915+105 microquasar could 
be explained in the framework of the EXORM (Stuchlik ct al., 2007d; Slany and 
Stuchlik, 2007). Therefore, we briefly summarize these results, and then we con- 
sider the case of Sgr A* , assuming relevance of all the three frequencies reported by 
Aschenbach (Aschenbach, 2004, 2007). 

4.1 GRS 1915+105 

Assuming mass M = 14.8 M and dimensionless spin a = 0.9998 for the GRS 1915+105 
Kerr black hole (see Stuchlik et al., 2007d; Slany and Stuchlik, 2007 for details), 
the EXORM predicts the following pattern of observable frequencies composed from 
the humpy and epicyclic frequencies and their combinations 



The corresponding humpy radius is = 1.29. At such a radius, the vertical 
epicyclic frequency of a particle orbiting the Kerr black hole with the mass and 
spin given above reaches the value v v = (0.23 ± 0.01) kHz that enters the highest 
(combinational) frequency. 



v x ~ (v T - u h ) = (26 ± 2) Hz , 

u h = v 2 = (41 ± 1) Hz , 

v T = u 3 = (67+ l)Hz, 

f 4 - {vi + Vh) = (108 ± 2) Hz , 

v b - (i/ v - v T ) = (0.17 ±0.01) kHz. 



(24) 
(25) 
(26) 
(27) 
(28) 
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4.2 SgrA* 

There are three frequencies related to the SgrA* QPOs (Aschenbach, 2004; T6r6k, 
2005): 

v u = 1.445 niHz , v m = 0.885 mHz , vy = 0.459 mHz . (29) 

These frequencies come in the rational ratio is u : v m : v\ <~ 3:2:1. Assuming EXORM 
with the humpy and radial epicyclic frequency and their combinational frequencies, 
we can distinguish three different cases of the resonant phenomena explaining the 
observed frequencies. 



v T : ~ 3:2 

The observed frequency pattern is given by u u = u T , v m = i/^, v\ = v T — iv The 
value of the spin, given by the frequency ratio, and the humpy frequency are then 
given by (compare Figs 2 and 3) 

a 3 :2 = 0.999984 , u h(3:2) = 590 ^ Hz . (30) 

Using the condition v m = i/h, we obtain mass of the black hole to be 

M = 0.667 x 10 6 M Q . (31) 

(Note that in this case we have chosen - see Fig. 3 - the value of the spin at the 
region where the frequency ratio v T : starts to be close to the asymptotical value 
of- 3:2.) 



u r : v\y ~ 2:1 

The observed frequency pattern is given by v u = v T + i/h, v m = v T , v\ = Vh- The 
value of the spin and the humpy frequency are then given by (compare Figs 2 and 3) 

ci2:i = 0.99925, i^ ( 2 : i) = 567 Hz . (32) 

Using the condition v\ = u^, we obtain mass of the black hole to be 

M = 1.235 x 10 6 Mq . (33) 



v T : fh ~ 3:1 

The observed frequency pattern is given by ^ u = ^ r ,^ m = v x — v^v\ = v^. The 
value of the spin and the humpy frequency are given by (Figs 2 and 3) 

£i3 : i = 0.999984 , i/ h(3:1) = 425 ^ Hz . (34) 



464 Z. Stuchlik, P. Slany and G. Tdrdk 



Using the condition v\ = v^, we obtain mass of the black hole to be 

M = 0.926 x 10 6 M . (35) 

However, by analysing data of the orbits of stars moving within 1000 light hours of 
the Sgr A* black hole, its mass is estimated to be M ~ 3.6 x 10 6 M©, and the error 
of the estimate is given by (Ghez et al., 2005; Weinberg et al., 2005) 

2.8 x 10 6 M Q < M < 4.6 x 10 6 M . (36) 

Clearly, the Sgr A* black hole mass predicted by the three simple variants of the 
EXORM are all completely out of the mass estimates given by relatively precise 
star-orbit measurements. Therefore, we can conclude that the EXORM cannot be 
applied to explain the QPOs observed in the Sgr A* black hole candidate. Although 
we expect a fast rotating black hole in the Galaxy Centre, it is probably not a near- 
extreme Kerr black hole that sites in Sgr A* . 



5 CONCLUDING REMARKS 

The extended orbital resonance model with the hypothetical humpy oscillations 
could be related to the QPO resonant phenomena in both thin Keplerian discs 
and marginally stable tori orbiting near-extreme Kerr black holes. The non-linear 
resonance is assumed between oscillations with the humpy frequency and the radial 
or vertical epicyclic frequency. Both parametric and forced resonance phenomena 
are possible, therefore, the combinational frequencies are allowed too. Generally, 
more than two observable oscillations are predicted. 

The EXORM can successfully explain all five QPO frequencies observed in the 
microquasar GRS 1915+105 (Stuchlik et al., 2007d), where a near-extreme black 
hole with a <~ 1 is predicted by spectral X-ray continuum fitting (McClintock et al., 
2006). Although Middleton et al. (2006) refer to a substantially lower, intermediate 
value of the black-hole spin, a ~ 0.7, to which the EXORM cannot be applied, the 
recent analysis by Narayan et al. (2007) and McClintock et al. (2007) demonstrated 
convincingly that the near-extreme black hole is more probable there, making the 
predictions of the EXORM still well viable. This model is successfully applied to 
explain QPOs observed also in the case of the other microquasar XTE J1650— 500, 
and in the near-extreme intermediate-mass Kerr black hole candidate in the system 
NGC 5408 X-l (Slany and Stuchlik, 2007). 

On the other hand, the EXORM is not able to explain, in any of its variants, 
the QPOs observed in the Sgr A* source, i.e., in the Galactic Centre Kerr black 
hole. The differences in the black hole mass estimated by the EXORM and by 
analysing the star orbits in vicinity of the black hole are too high to believe that 
some more exact QPOs measurements could give a solution of this discrepancy. On 
the other hand, it is worth to note that the Sgr A* QPOs, if true, could be explained 
within the multiresonance model based on the assumption of strong resonances of 
oscillation with the orbital and both epicyclic frequencies (Stuchlik ct al., 2007a). 
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The spin estimated in this case is <~ 0.983, corresponding to a fast rotating, but not 
near-extreme Kerr black hole. 

We can conclude that the EXORM could be considered as a promising model 
of QPOs in near-extreme Kerr black holes systems, where the oscillations occur in 
the innermost parts of the accretion disc. The model enables very precise measure- 
ments of the black hole parameters, in particular, of the black hole spin since its 
value is given (usually) by the frequency ratio of the humpy frequency and the radial 
epicyclic frequency. However, the predictions of the EXORM have to be confronted 
to the black hole parameter estimates coming from other methods, as those based 
on the optical phenomena in strong gravitational fields, e.g., the spectral X-ray 
continuum fitting, the profiled spectral (Fe-K) lines fitting, time delay methods, or 
by methods based on the orbital motion analysis. The present situation seems to 
be rather controversial, however, we believe that in future, with developing both 
the theoretical models and observational techniques, we could be able to under- 
stand the accretion phenomena in much deeper detail as compared to the present 
understanding. 
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